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*y-\ ' Abstract: We prove the theorem announced in Phys. Rev. Lett. 85:5022, 2001 

^^ \ concerning the existence and properties of fractal states for the Schrodinger 

[^. ■ equation in the infinite one-dimensional well. 

o 
o 



1. Introduction 

Fractals are sets and measures of non- integer dimension p9[ . They are good 
models of phenomena and objects in various areas of science. They are often 
connected with non-equilibrium problems of growth [|22| and transport |llj|4|. In 
this direction their importance can be expected to grow in view of the recent re- 
sults connecting transport coefficients with fractal properties of hydrodynamic 



modes [[13,12 . Their ubiquity in dynamical systems theory as attractors, re- 
pellers, and attractor boundaries is well-known P3|,p|. 

Fractals have also been found in quantum mechanics. For instance, quan- 
tum models related to the problem of chaotic scattering often reveal the fractal 
structure p,[jj],f|] .Quantum field theories in fractal spacetimes have recently 
been reviewed in fllcf . 

Recently fractals have been shown to satisfy Schrodinger equations for sim- 
plest non-chaotic potentials 0,pq|. In this paper we prove some of the results 
announced in pqj . 

This work has two objectives. One is to present the first (as far as we know) 
rigorous proof of fractality of any quantum states. The other is to illustrate a 
convenient method of calculating dimensions of graphs of continuous functions 
introduced by Claude Tricot p4|. 
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One of the oldest fractals is Weierstrass function p5|,p|: 

oo 

W{x)=^a n cos{b n xTi). (1) 

introduced as an example of everywhere continuous nowhere differentiable func- 
tion by Karl Weierstrass around 1872. Maximum range of parameters for which 
the above series has required properties was found by Godfrey Harold Hardy in 
1916 urn, who also showed that 

sup{\f(x)-f(y)\:\x-y\<S}~5 H , (2) 

where 

ln(l/a) 

From this it easily follows (see below) that the box-counting dimension of the 
graph of the Weierstrass function W(x) is 



mo 



In a 



In 6 



(3) 



Functions whose graphs have non-integer box-counting dimension will be called 
fractal functions. Even though the box-counting dimension of the Weierstrass 
function is easy to calculate |2J], the proof that its Hausdorff dimension has 
the same value is still lacking. Lower bounds on the Hausdorff dimension of the 
graph were found by Mauldin [glJgCJ] • Graphs of random Weierstrass functions 
were shown to have the same Hausdorff and box-counting dimensions for almost 
every distribution of phases |lq| . 

The above construction can easily be realized in quantum mechanics. Consider 
solutions of the Schrodinger equation 

id t V{x,t) = -V 2 <f(x,£) (4) 

for a particle in one-dimensional infinite potential well. The general solutions 
satisfying the boundary conditions ^(O, t) = = \P(tv, t) have the form 

oo 

iP(x, t) = ^ a n sin(nx)e~ m2 \ (5) 

n=l 

where 

2 r 

a„ = — / dx sin(nx)^ r (x, 0). (6) 

n Jo 

Weierstrass quantum fractals are wave functions of the form 

M 

& M (x, t) - N M Y, <? n(8 ~ 2) sin^aOe-^"', (7) 



n=0 



where q = 2, 3, . . ., s€ (0,2). 
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In the physically interesting case of finite M the wave function \&m is a solution 
of the Schrodinger equation. The limiting case 

oo 

W(x,t):= Van V M {x 1 t)=N S jT q n{ - s -^sm(q n x)e- lq2n \ (8) 

n=0 



with the normalization constant N = \ — (1 — g 2 ( s 2 )), is continuous but nowhere 

differentiable. It is a weak solution of the Schrodinger equation. Note that (ph 
converges for |g i_2 | < 1 = s < 2. Later on we will show that the probability 
density of (g) shows fractal features for s > 0. Thus the interesting range of s is 
(0,2). 

We show that not only the real part of the wave function \&{x, t), but also 
the physically important probability density P(x,t) := |i^(2;,i)| 2 exhibit fractal 
nature. This is not obvious, because ^(a;, £)| 2 is the sum of squares of real 
and imaginary part having usually equal dimensions. One can easily show that 
the dimension of the graph of a sum of functions whose graph have the same 
dimensions D can be anythingn from 1 to D. 

Our main results are given by 

Theorem 1. Let P(x,t) denote the probability density of a Weierstrass-like wave 
function m). Then 

1. at the initial time t = the probability density, Pq(x) = P(x,0), forms a 
fractal graph in the space variable (i.e. a space fractal,) of dimension D x = 
max{s, 1}; 

2. the dimension D x of graph of Pt(x) = P(x,t = const) does not change in 
time; 

3. for almost every x inside the well the probability density, P x (t) = P{x = 
const, t), forms a fractal graph in the time variable (i.e. a time fractal,) of 
dimension D t (x) — D t := 1 + s/2; 

4-. for a discrete, dense set of points Xd, Px d (t) = P(Xd,t) is smooth and thus 
Dt(x d ) = 1; 

5. for even q the average velocity —&{t) is fractal with the dimension of its 
graph equal to D v — max{(l + s)/2, 1}; 

6. the surface P{x,t) has dimension D xy = 2 + s/2. 

The physical meaningof the above theorem has been discussed in J26| . We 
give its proof in section |3|. It is based on the power-law behavior of the aver- 
age (S-oscillation of the infinite double sum present in P(x,t) = ^(a;,^)! 2 (|55|). 
In section @ we introduce some fundamental concepts and facts to be used in 
subsequent analysis. Calculations of probability density and average velocity are 
provided in the Appendix. Positive real constants are denoted by c, c±, C2, . . . 



2. Box-counting dimension 

In this section we recall several equivalent definitions of box-counting dimension, 
state a criterion for finding the dimension of a continuous function of one variable 

1 Let /i and J2 be functions with graphs having dimensions, respectively, 1 < D\ < D2 < 2. 
Let gi = /1 + /2, 92 = /1 — /2- Then the box-counting dimension of both the graph of g\ and 
gi is -D2, but the dimension of the graph of their sum 91 + 92 = 2/i is D\ 6 [1, Z3a[. 
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and prove a connection between the dimension of graph of function of n variables 
and the dimensions of its sections. All of this is known with the exception, 
perhaps, of theorem |3j which might be new. We concentrate on the theory of 
box-counting dimension for graphs of continuous functions of one variable. More 
general theory and a deeper presentation can be found for instance in [9, 10, 24fl. 



Let A C R™ be bounded. Consider a grid of n-dimensional boxes of side 5 

[miS, (mi + 1)5] x ... x [m n 5, (m„ + 1)5]. (9) 

Let N(6) be the number of these boxes covering the set A. It is always finite 
because A is bounded. 

Definition 1. Box-counting dimension of the set A is the limit 

dim B (A): = lim^. (10) 

If the limit does not exist one considers upper and lower box-counting dimensions 

-i — / ,s , lnN(5) , , 

dim B (A) := hmsup— - — — , (11) 

s^q In 1/6 

dim^^^hminf^T? ( 12 ) 

5 — ►o In l/o 

which always exist and satisfy 

dim R (A) > dim B (A). (13) 

The box-counting dimension exists if the upper and lower box-counting dimen- 
sions are equal. 

Several equivalent definitions are in use (see |IJ,y,|lO,24| for a review). The 
most convenient definition to study the fractal properties of graphs of continu- 
ous functions is given in terms of 5- variations p4[ . It is essentially a variant of 
Bouligand definition [El . We shall restrict our attention to dimensions of curves 
being subsets of a plane. 

Let Ks(x) be a closed ball {y £ M 2 : \x — y\ < 5}. 

Definition 2. Minkowski sausage or ^-parallel body of A cM. 2 is 

A s := |J K s (x) (14) 

xeA 
= {y£U 2 : 3x£ A, \x-y\<5}. (15) 

Thus the Minkowski sausage of A is the set of all the points located within 5 
of A. 

Proposition 1. The box-counting dimension of a set A CM 2 satisfies 

*-<^]5s( 2 -^). ™ 

where V(5) = vol (As) is the area of the Minkowski sausage of A. 
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Proof. Every square from the 5-grid containing x € A is included in Kr^ s {x). 
On the other hand, every closed ball of radius v2<5 can be covered by at most 
16 squares from the grid. Therefore 

5 2 N(e) < V(A^ S ) < m 2 N{e). (17) 

□ 

Consider a continuous function on a closed interval / : [a, b] — > M. Its graph 
is a curve in the plane. To find its box-counting dimension estimate the number 
of boxes N(6) intersecting the graph. 

Choose column {(x,y) : x € [nS, (n+ 1)5]}. Since the curve is continuous the 
number of the boxes in this column intersecting the graph of / is at least 



sup f(x) — inf f{x) 



/S (18) 



and no more than the same plus 2. If / was a record of a signal then the difference 
between the maximum and minimum value of / on the given interval quantifies 
how the signal oscillates on this interval. That's why it is called oscillation. 

Definition 3. 5-oscillation of f at x is 

osc s (x)(f) := sup f(y)- inf f(y) (19) 

\y-x\<6 \y-x\<S 

= sup{|/(y)-/(z)| : y,z€[a,b]n[x-d,x + 8}}. (20) 

We will skip (/) if it is clear from the context which function we consider. 

From ( |18| ) we obtain the following estimate on the total number of boxes 
covering the graph of /: 

M M 

J2 osc 5/2 (x m )/S < N(S) < 2M + J2 osc s/2 {x m )/S, (21) 

m— 1 m— 1 



where x m = a + (m — 1/2)6 is the middle of the m-th column from the cover of 
the graph and M = [^r^l is the number of columns in the cover (|"af| stands for 
the smallest integer greater or equal to x). Thus 

N(S) w Moscg/2/S. (22) 

If the graph of / has the box-counting dimension D, N(S) scales as S~ D . This 
implies the following scaling of the oscillations 

OSC5/2 ~ N(S)S/M ex S 2 - D . (23) 

We have thus suggested a connection between the box-counting dimension of the 
graph and the scaling exponent of the average oscillation of the function /. 

Definition 4. 6- variation of function f is 

Va r<5 (/):= f osc s {x)(f)dx (24) 

J a 

=:(b-a)osc 5 (f). (25) 
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Geometrically, variation is the area of the set scanned by the graph of / moved 
horizontally ±o~ and truncated at x — a and x = o, thus it is a kind of Minkowski 
sausage constructed with horizontal intervals of length 26. This observation leads 
to a convenient technique for calculating dimensions. 

Theorem 2. Let f(x) be a non constant continuous function on [a,b], then 

dim* graph/ = lim f 2 - ^^V (26) 

(5^o \ mo J 

The proof consists of showing equivalence of Var^/) with the Minkowski sausage 
and follows from inequality (Q, p. 130-132, 148-149) 

Var,(/) < V(A 5 ) < cVar,(/), (27) 

where 

A = graph/ 
C= Cl + c 2 /s 

sup f(x) - inf f(x) 

c£[a,b] xe[a,b] 

This is where the assumption of non-constancy of / comes in. Derivation of ( p7J ) 
is not difficult but rather lengthy and will be omitted. 

This theorem is the main tool to prove Theorem |]. In order to find the 
dimensions we will look for estimates of ^-variation. They will usually take the 
following form: 

Proposition 2. 

1. oscs{x)f(x) < c5 2 ~ s => dims graph/ < s. 

2. W := J* \f(x + 5)- f(x - S)\dx > c5 2 - s => dims graph/ > s. 

Proof. 

1. Vara/ = J* osc s {x)(f)dx <{b- a)c5 2 - s . 

2. osc 2& {x)f > \f(x + S)- f{x - S)\ => Var 5 / > (b - a)c{5/2f- s . 

D 

To prove the last point of theorem fl\ we need to know what is the dimension of 
the graph of / : R" — > M given all the dimensions of its one-variable restrictions. 

Theorem 3. Let f € C°([ai, &i]x. . .X [a„, 6„]). For every point x — (x 1 , . . . ,x") g 
[ai, fei] x . . . x [a n , b n ] define x 1 :— (x , ■ . ■ , x 1 ^ 1 , x 1+1 , . . . , #"). TTiera 

fi [x \ (x ) := f (x Q , . . . , Xq ,x,Xq , . . . , x ) (28) 

is a restriction of f to a line parallel to i-th axis going through xo and /i[£o] *= 
C°([oi,6i]). 

i. //Vx : oscsfiix 1 } < CiS Hi then dims graph./ (a; 1 , . . . , x n ) < n+l—mm{H\, . . . , H n }. 
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2. IfVax$fi\xQ\ > Ci6 Hi for a dense set Xq € A C A = [ci, 61] x . . . x [a,;_i, bi-%] x 
[a l+ i,6 4+ i]x. . .x[a n ,6„] £/ien dims graph/ (a 1 , . . . ,x n ) > n+l-mm{Hi, . . . ,H n }. 

3. If all of the above conditions are satisfied then 

dims graph/ (a; , . . . , x") = n + 1 — min{i?i, . . . , i/„} 

= n — 1 + maxjsi, . . . , s n } 

where 

Si = sup dims graph. fi[x l ](x l ). 

X 1 

In other words, the strongest oscillations along any direction determine the 
box-counting dimension of the whole n + 1-dimensional graph. 

Proof. We will show the theorem for n = 2 for notational simplicity. General- 
ization to arbitrary n is immediate. Let / : [ai,&i] x [02,62] —> K. Divide the 
domain into squares Xi x Yj of side 8. This gives rise to K columns Aij of <5-grid 



(61 — ai)(&2— 0-2) 

1 . The number of 5-cubes having a common point with the graph of / in column 
A.^ is not greater than (sup^. . / — inf ^ . f)/S + 2. But 

\f(xi,yi) - f(x2,V2)\ = \f(xi,yi) - f(xi,y 2 ) + f(xi,y 2 ) - f(x2,V2)\ 

< \f{xi,Vi) ~ /(zi, 2/2)1 + \f{xi,y 2 ) - f(x 2 ,y 2 )\. 

Therefore 

sup/-inf/ = sup \f(xi,yi)- f(x 2 ,y 2 )\ 

< sup sup f(x,y) + sup sup f(x,y) 

igXi ySYj yGYj x£Xi 

< sup osc s/2 fi[x] + sup osc s /2 f 2 [y] 

x£X t yeYj 



Thus 



< c< 5 min { H i^ 2 }. 

\n{Kc8 mi ^ H ^ H ^/8) 



dims graph/ (x , x ) < lim 



<5->o In 1/8 

<3-rxun{Hi,H a }. 

2. Set x S Xi. From (|l7|) and @) it follows that the number Ni(8) of 5-cubes 
in columns Aij covering the graph of f 2 [x] (y) and the variation of f 2 satisfy 

Vaisf2[x]<cS 2 Ni(6). (29) 

Thus 

Ni(6) > c sup Vnr s f 2 [x}/5 2 > cS H2 - 2 . (30) 

xeXi 

Therefore the number N(8) of boxes covering the whole graph of / satisfies 

M 

N(8) > cj] sup Vai s f 2 [x}/8 2 > Cl S H2 - 3 . (31) 



i=l 
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The same can be repeated for any direction thus 

N(8) > C2 5 min ^ Ht ' H ^- 3 . (32) 

3. An immediate corollary. 

Generalization to arbitrary n is achieved by observing that KS n ~ const. D 

Remark 1. Another definition which has some convenient technical properties is 
the Hausdorff dimension fT^,R,p), however it is often too difficult to calculate. 
For instance, as mentioned before, there is still no proof that the Hausdorff 
dimension of the Weierstrass function is equal to its box-counting dimension. 
Thus in practice one usually uses the (upper) box-counting dimension. This is 
also our present approach. It is often assumed that the box-counting dimension 
and the Hausdorff dimension are equal. A general characterization of situations 
when this conjecture really holds is also lacking. 

3. Proof of Theorem @ 

1., 2. We will show that for every fixed t the graph of the probability density 
\\P\ 2 (pq) as a function of x has the box-counting dimension s. 
(a) Fix t. Let 

n k 

P n ( x ) := J2 Q k{s ~ 2) J2 sinq l xsmq k ^ l x COs(q 21 - q 2( - k -^)t, 

k=0 1=0 

q = 2, 3, .... It is a smooth function whose derivative at every point satisfies 

n k 



k ~ l x\ 



where 



\K{x)\ ^^(^^Icos^sinq 

fc=0 1=0 

n fc+1 

< 2Yq k(s - 2 ^ < d 1 (s ) q)q^ 8 -^ 

^-^ (7—1 

k=0 y 

di(s,q) 



(q-l)(qs-i-i)- 
Let S = q- n . Then 

osc s (x)P n < 26 sup \P^(x)\ <2d 1 {s 1 q)5 2 ~ s . (33) 

se[0,7r] 

On the other hand, for 

R n {x) ■- P{x) - P n {x) 

oo k 

= V" q k (s-2) Y sing 1 x sin q k ~ l xcos(q 21 - q 2< - k ^)t 

k=n+l 1=0 
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we have 

osc 5 ( X )R n <2 £ g «- a >(fc + l)< 7 " . 

k=n+l ^ q ' 

Polynomial growth is slower than exponential, therefore for arbitrarily 
small e there is some M such that 

Vn > M : n< (q £ ) n . 

This leads to the following estimate of the oscillation of R n : 

osc s (x)R n < d 2 (s,q)5 2 - s - e , 

where 

<fe(s,g) 



iq 



8—2 



(l_g«-2)2- 

Thus for all x and 8 = q~ n , where (lnn)/n < ehiq, we have 

osc$(x)P < oscs(x)P n + oscs(x)R n 
< (2d 1 +d 2 )6 2 - s ~ £ . 

From proposition it follows that 

dims graph P t (x) < 2 — (2 — s — e) = s + e — ► s. 

(b) Fix t. Let f(x) — P(x,t). We want to show that 

,6 

W := / \f(x + 5)- f{x - S)\dx > cS 2 - s . 

J a 

Take a — 0, b — tt. Notice that (we skip the normalization constant) 

W= f : \f(x + S)-f(x-S)\dx (34) 



J2 q k( - s -V J2{smq l (x + 6) smq k - l (x + S) + 
k=0 1=0 

sinq l (x — S)sinq k ~ l (x — 6)}a,ki\ dx 

oo k 

^ g fc(s-2) "^{cos q l xs [ n q k - l xs - m q l $ cos q k ~ l 6} a k i dx, 



k=0 



1=0 



where 



a ki = cos(q zk - q M )t 
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Take \h(x)\ < 1. Observe, that 



£*(*) 



dx > / |/i(a;)| 



£/*(* 



rf.i: 



> 



> 



y^h(x)fj(x)da 



(35) 



/ h(x)f k (x)dx 

J a 


-£ 


/ h(x)fi(x)dx 

J a 



One can interchange the order of summation and integration because f{x) 
is absolutely convergent. Let us take S — q~ N , h(x) — smq m xcosq n x. 
After substitution in (p3) using (F55|) we obtain 



W > 



£) 9 fc(-2)£; sin ^ sing 



k-l-N 



fc=0 



i=0 
J_„ k-l 



dx sin g x cos q ' x cos q m x cos g n x a^ 



u 



= _„("*+«)(s-2) 



}-AT 



-iV 



cos q sin g cos 



3 ( g 2(m+n) _ g 2m)j 



In what follows we will show three consecutive proofs that 



3c: W = g (m+n)(s_2) 



cos q m N sin 



ig"- Ar cos(g 2(m+ " ) -g 2m )i 



>cq N ^ s ~ 2 \ 



for £ = ^f , t/n £ Q, and a general proof for arbitrary real t. We will take 
advantage of the fact that q is integer. 
Let TV = m + n. Then 

W = q N{s - 2) \cosq m - N smq- m COs(q 2N - q 2m )t\ . 

It is enough to consider t £ [0, ir[. 
i. Let t — ^-. Take m such that 2m > I, for instance ra — l. Then 



7 



(^.^^(^-J.^ 



therefore 



cos 



q l 



= 1. 



We also have sin q l = const and 



l-N 



cos 1 = cos q < cos q < cos g 



cosO. 



This means that for t = ^f , for sufficiently large N 



W > q (" 2 ' cos 1 • sing • cosO = const g 



N(s-2) 
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ii. Let t/n G Q. For the next two proofs let us write t/Tt in q basis 





t 

IT 


fll 

q 


+ 


«2 a 3 


fc=i 


q k 


re a k <E {0, 1, 


..., 


Q- 


1} 


Therefore 






cos[(q 2N - q 2 


m )A 












= cos^g 2 ^ 


"V 


+ (> 


f 2JV 


~ 2 a 2 + ...+ 


02w + q 


-1 



(36) 



-<? 2m -v 



02W+1 + ■ 

7 2m ~ 2 a 2 + . . . + a 2m + q~ 1 a 2m +i + ■■■)} 



0-2N+1 — 02m+l Q-2AT+2 _ 02 m +2 



9 



9 



(37) 



If we could only choose to so that the first two terms in this series 
cancel out, we would have a lower estimate on the cosine, because in 
this case 



02AT+3 — 02m+3 



Thus 



cos[(g 



2 A 



-"'-"i^ 



q 2m )t] > cos(tt/o 2 ) > costt/4. 



1 



Take a rational t. In any basis q the expansion (|36|) of t/n is finite or 
periodic. Finite expansions have been treated in the previous point, 
thus we assume here the expansion of t/ir is periodic with period T. 
Thus for sufficiently large n we have 



On+T — o„ 

so that t/ir can be written as 

t/n = 0.ciia 2 . . . a K (a K+ i . . . a K+T ), 

where (a«-+i ■ ■ ■ clk+t) means "repeat cik+i ■ ■ ■ clk+t periodically ad 
infinitum" . We shall now estimate W for this t. Every N > K can be 
written as N = K+pT+r, wherep is natural or and re {1,2,..., T}. 
If we now take m = K + r then not only the first two but all the terms 
in (|37|) cancel out, therefore 



COs[(q 



2 A 



,2m 



)*] 



1. 



Obviously, 



Therefore 



sin<r m >sinq-^ +T \ 

~ N = cosq~ pT > cosl, 



cosq 



JV(s-2) 



(K+T) r 



W > q" (S ~ Z > sin q- {JX+1 > cos 1 = const q 



JV(s-2) 
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iii. General case of an arbitrary t/it G [0, 1]. Its expansion in q basis is 
given by ([36|) . Let A be the set of all the two element sequences with 
elements from the set {0, 1, . . . , q — 1}. Thus 

A = {{0,0},{0,l},...,{0,q-l},{l,0},...,{q-l,q-l}}, 

we write Ak,i := {k, I}, k, I € {0, 1, . . . , q — 1}. Consider all the pairs 
of consecutive g-digits of t/n of the form 

{02m+l, »2m+2}, (38) 

i.e. {ai, 02}, {123, 04}, etc. Every such pair is equal to some Ak,i- Let 
Nk,i be the first such m for which 

Ak,l = {«2m+l) 0>2m+2J- 

If Ajy for given fc, Z doesn't appear in the sequence of all the pairs ( p8| ) 
we set Nk,i = 0. Let 

M = supJV fc ,i. 

k,l 

Thus if n > M the sequence {a2 n +i, fl2n+2} has appeared at least once 
among the pairs {ai, a 2 }, {a 3 , a 4 }, . . . , {a 2 M+i, «2M+2}- Therefore, for 
every N > M we can find such an m E 1,2, ... , M that 



7T 



7T 7T V2 
> COS — > COS — = . 

q 2 ~ 4 2 



Also 



which leads to 



sin<T m >sin<T M , 
cosq m ~ N > cosq M ~ N > cosl, 



W > g^*- 2 ) ^E sin g- M cosl = const q N ^- 2 K 

We have thus shown that for every t, for natural q and for 6 = q~ N 

W > const 5 2 - s , 
therefore (proposition ||) 

dims graph Pt(x) > 2 — (2 — s) = s. 
3. For almost every x, D t (x) — dims graphP x (t) = D t := 1 + s/2. 



We will use the form (pq) of the probability density. It is enough to analyze 
the dimension of 

oo c 

P(t) := ]T 9 2c (- 2 ) sin^^g-^ 8 - 2 ) sing c - d zcos[(<2 2c - g 2 ^)*]. 
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(a) Let 



P n (t) :=J2<1 2C{S ~ 2) smq c xJ2 c r d(S ' 2) smq c - d xcos[(q 2c - q 2(c - d >)t}. 



Then 



\K(t)\ 



J2 1 2C(S ~ 2) sin q c x J2 q- d{s - 2) sin q c - d 3 



d=l 



.tfc _ q 2(c-d) )s[n[ ( q 2c _ q 2(c-d) )t] 



< j^ q 2c{s - 2) J2 q- d(s ~ 2) (q 2c - q 2(c - d) ) 



d=l 



^ g 2c( S -2+l)^ g -d( s -2) (1 _ g - 



2d \ 



J2-s 



q 2 ~ s - 1 



q~ s - 1 



d=l 
ns _ -i 

q s — 1 



2n(s-l) _ 1 



q 2(s-l) _ 1 



, g n(s - 2) -i 2 r s -i) g 2 " (s " 1) -i 



q s-2 _ I 



,2(8-1) _ 1 



Therefore, for n large enough 



\P' n {t)\ < dg" 111 ^ 8 ' 2 ^- 1 )' 8 - 2 } = eg" 5 . 
Let <5 = <T Q ™. Then q n = (T 1 /" and 

osc s (t)P n < Cl 26q ns = 2ck5 1 - s/q . 



Let 
Then 



i?„(i) :=P(t)-P„(t). 



OSC,(i)P„<2 J] q 2 ^- 2 ^^^ 

c=n 

2q 2 



c—n+1 d—1 

,2-s x 



< 



tE« 



2c(s-2)+c(2-s) 



c— n+1 
^(s-2)/a_ 



To obtain consistent estimate we must set 

s _ 2 s 
a a a 

which gives a = 2. Thus 

osc,5(i)P< (2c 1 +c 2 )(5 1 - s / 2 . 
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(b) Now we want to show that 



W 



dt\P{t + 6)-P(t-6)\ >c5 x - sl2 . 



Set a = 0, b = 2tt for convenience. Then 

f27T 



!! = / dt 
/o 



V^ g 2c(s-2) 



S inq c xJ2<l~ d{S ~ 2) 



smq c d x ■ 



c=l 



d=i 



{cos[(g 2c - q 2 ^- d ^){t + 8)] - cos[(q 2c - q 2 ^-^)(t - S)]} 

c 

y^ q -d( s -2) 



dt 



V\ 2c(s 2 ^sing"x 



c—d 

sing x ■ 



d=l 



• {-2sin[(q 2c - q 2{ - c "^)t] sin[(q 2c - q 2{c - d} )6]} 

Using our standard arguments we multiply the integrand by a suitable 
function smaller or equal to 1: 

/•27T _ 

W> / dt\h(t)\\P(t + 6)-P(t-5)\ 
Jo 

r 2ir 

• / dth(t)[P(t + 5)-P{t-5)} 



sm(q c x) sm(q c ~ d x) sin[(q 2c - q 2i - c -^)q- 2N ] 





We choose h(t) = sin[(<? 2c - q 2 ^ c - d ^)t\ and set 6 = q~ 2N . It follows that 
W > 2^ 2c - d ^ {s - 2 ^ 

We now want to show that for almost all x 

W > CzS 1 -^ 2 = c 3 q N ( s -V- 
Set 2c - d = iV. Then 

W > 2 7 rg JV(s - 2) sin(g c a;) sin^^x) sinfe 2 ^-^ - q~ 2c 
Thus it is enough to bound 

/ N — r \ ■ r 9fV_N"l 



sm(q c x)sm{q N - c x)sm{q 2l - c - N '> ■■ J ' ' 



(39) 



(40) 



from below. 

Choose rational x/tt. All the rational numbers in a given basis q have one 
of the two possible forms: finite or periodic. In the first case (x/tt = k/q l ) 
we cannot find the lower bound on (|o|). We cannot succeed, because at 
these points the function P x (t) is smooth (cf . the next point of the proof) . 
The other case means that x/tt can be written as 



x/tt = 0.aia 2 ■ . ■ a K (a K +i ■ ■ ■ a K+T ), 
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where again (clk+i ■ ■ ■ ik+t) denotes the periodic part. Therefore, for ev- 
ery N > K, q n xmodir can take only one of T values: q K+1 xmodir, . . . , q K+T xmodTT. 
Let us take c = 1, N > K . Then | sm(<7 c a;)| = | sin(gx)| > and is a con- 
stant. | sin(g Ar_1 x)| takes one of T values, none of which is 0, therefore it 
is always bounded from below by 

inf \smiq K+l x)\ > 0. 

1=1,2,.. .,T 

Also the last term can be bounded: 

| sin(cT 2(Ar - c) - q- 2c )\ > sm(q- 2 - g" 2 ^" 1 )) > sin^ 3 
for N > 3. Thus for rational x/tt with periodic expansion in q 

W>c 3 q N (°- 2 \ 

where C3 = 2tt\ sin(gx) sin(g _3 a;)| mfi=i,2,... t | s'm(q K+l x)\. 
Consider now irrational x/tt. Inequality (|3^) for c = N takes form 

W > 2irq N ( s - 2) |sin(g Jv x)sinxsin[l - q- 2N }\ 

>cq N ^-^\sm(q N x)\, (41) 

for N > 2. Instead of showing it can be bounded from below we will use 
it to prove that for almost every x 

dim B graphic (i) > 1 + a/2. (42) 

(43) 



Let 


x n '■= Q n (x/7r) mod 1. 




Let 






Fn ■= {* 


:3n>N [x n < -^j V (l - 


Xn - q Na 


a e [0,1]. Let 


OO 

*s ■= n f n 





Clearly, 

*N -> *N+1 -> *N+2 ■ ■ ■ 

Since Renyi map ( |43| ) preserves the Lebesgue measure we have 

til%) < 2 ( -4- + -j^- + ~^ + ■■)= -^i ( 44 ) 

\q"a q" +l a q" +z a J q — lq" a 

Therefore 

0< M (F£)<mf M (F£) = (45) 

It follows that for almost every {x n } 



-na 



In sin x In o n 
lim — > lim > lim > — ahxq. (46) 

n — >oq Tl n — »oo 77, n — >-oo 2?2 
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Thus for every a > and for almost every x/tt G [0, 1] we have 

dim B graphic) = lim (2 - ^p^ ) (47) 

>2 + limf lnVa ^(^ (48) 

But 

VmPx(t) > W (49) 

therefore from (Edl) 

, „ / n „ /In c + iV(s — 2) lno + In I sin(a;„7r)|\ , „. 

dim B graphP x (t > 2 + lim ^ ^— 2 ! ^-^ 50 

\ 2N mq J 

/In I sin(x„7r)| \ 
= l + s/2 + llm ( I^JI j (51) 

> 1 + ,s/2 - a. (52) 

But a is arbitrary, thus 

dim B graphP x (£) > 1 + a/2. (53) 

4. For a discrete, dense set of points Xd, D t {xd) — dims graphP^ (t) = 1. 

Let Xk, m — t 11 , where k G N, m = 0, 1, . . . , q k — 1. The set {£fc ]TO } is dense 
in [0, 1]. At these points ^(xfc, m ,t) is a sum of a finite number of terms: 

'fe-i 



* Of-V = V ' (' " F^O & ( -**r'-i^ 



Therefore, 

,. , ,- , m-K 

dims graph 






= 1. 



5. For even q i/ie average velocity -gj (t) is fractal with the dimension of its 
graph equal to D v = max{(l + a)/2, 1}. 

Heuristically, this is rather obvious because 

f(*)«EV sin * 2fc ' - E<^ 3)/2 ^fl**- 
fe=i y fe=i 

Thus the average velocity is essentially a Weierstrass-like function and the 
dimension of its graph should be 



It is enough to consider 



oo k(s—l) 

w r (t):=Efsrn 8fal fe a *- 1 >*- 

fc=l y 
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(a) Let 



W, 



JL gfc(*-i) 



'-(*)== ES*ri 8in (« afc - 1 )*- 



fe=i 



Set 5 = g- an . Then 



|W^(t)| 



^(Z^-^cos^-l)* 



fe=i 



<^ g fc(-l) < Cl< j(l-)/c 



fc=l 



where 

Therefore 
Now, for 
we have 



ci = 



9 s - 1 - 1 



os C(5 (£)VK„ < 2c l S < - 1 - s ^ a S = 2c l 5 l+il - s ^ a . 
P n (t):=W(t)-W n {t) 



\Pn{t)\ 



00 „fc(«-l) 



fe=n+l 



q 2k_ 1 



< 



X gfc(s-l) 



E 

fe=n+l 



~ fc(s-l) 

<2 ^ «-__ =C2 j- 



fc=n+l 



where 



Thus 



,2/,' 



f'2 



2g 8 



l-« 



s-3 ' 



oscs(t)P n < 2c 2 5 <* . 
To obtain consistent estimate for both P n and W n we must set 

1 — s 3 — s 



1 



o 



thus a = 2 and 6 = q . Therefore 



osc s (t)W < osc s (t)W n + osc s (t)P n < 2(ci + c 2 )<5 2 - V 



(b) Consider 



|W(t + £)- W(i- J)|di 

/b oo k(s— 1) 

dt E Vt cos (? 2fe - *)* s[ ^i 2k - l ) s 
fe=i q 



> 



dth(t)f N (t) 



E 



dth{t)f k {t) 
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where 



hin=^^ T cos(q 2k -l)tsm(q 2k - l)fi. 



q 2k_l 

Let h{t) = coa(q 2N - l)t, S = q~ 2N . Then 

rb 

h(t)f N (t)dt 



JV(s-l) 



-,2N 



-sin{l-q- 2N ) f cos 2 (q 2N -l)tdt 

^ J a 



>^_sin- / cos 2 (q 2N -\)tdt 
6 ' 



2iV 



q 

> l q ms-3) 



- 2 



b-a sin(2b(q 2 " - 1)) - sin(2a( g 2JV - 1)) 
~1T + A{q 2N - 1) 



2 



> I,jV(6_ ). 

"8 v ' 



b — a 


2-2 


2 4 
6- a " 


, g 2iV_ 


2 





On the other hand 

rb 



h(t)f k (t)dt 



Ms-1) 



,2fc 



- sin [((7 2fe - l)q- 2N ] / cos^ - l)icos(cz 2JV - l)t 



- l)tdi 



7 fc(«-i) 



< 



sin(6(^ iv - q 2k )) - sm(a(q™ - q 2k )) 



q 2k _ i 2 (q 2N - q 2k ) 

sm(b(q 2N + q 2k )) - sm(a(q 2N + q 2k )) 



2(q 2N +q 2k ) 



< 2q k(s - 3) 

< 2q k(s - Z) 

< 5 ■ q k{s - 3) 5 



\q 2N -q 2k \ q 2N + q 2k 



q 2N _ q 2(N-l) q 2N 



Therefore 



W > -S^ 3 - S ^ 2 (b -a)-J2 $q k{s - 3) S 



"I n s ~ 3 



>L S ( 3 -°)/ 2 (b-a)-5- 



,8-3 



-S. 
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But ^p < 1, thus for large enough N (small enough 8) the first term 



2 

dominates the other, therefore 



W>c6 2 - (1+s V\ 



with c = ^j^j for example. 
From theorem it follows that 



D v = 1±^. 
2 

6. The surface P(x,t) has dimension D xy = 2 + s/2. 

Setting iori constant we have shown that oscillations are bounded by c8 H 
where exponent H is one of l,s, s/2. We also showed the lower bound of 
variation is always c8 H again with H being one of 1, a, s/2. What's more, 
there is a dense set of points x for which V&TsP x (t) > c8 s ^ 2 . One can take for 
instance all rational x/tt with periodic g-expansion. Thus from theorem we 
have 

D xt = l+max{D x ,D t } = 2 + -. (54) 

D 



A. Appendix: Auxiliary calculations 

A.l. Probability density. Take the fractal wave function (pi) 



- oo 



f (s, t) = W - (1 - g 2 ( s - 2 )) ^ q"( s - 2 ) sin(g n x)e 
Let us calculate two useful forms of the probability density P(x, t) 

P(x,t) = \*(x,t)\ 2 



1 'l-q 2 ^- 2 ^ Y, q {m+nKs ' 2) smq n x S mq m xe~^ 2 "- c ' 2 ^ t . 

m.n—O 



Taking k = m + n, I = n we obtain 

k 



„-*(«" -V""")* 



P(x,t) = - (l - q 2{s - 2) ) Y^^ £) sing's sin g*-W 

k=Q 1=0 

oo fe 

- (l - q 2{s - 2) ) Y 1 k{S ~ 2) Y 8illqlxsinqk ~ lxCOS ( q21 ~ l 2 ^^)^ (55) 



fc=0 1=0 

Substitute c = m, d = m — n 
2 



l x+ 



P(x,t) = - (l-g 2 ^ 2 ') Y, {g 2m(s_2) sin 2 q r - 

771—0 

+2 ^ g (m+n)(«-2) s i ng « a;s i ng ™ a ; C os(g 2m - q 2n )t, 



n<m 
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- (l - q 2{s - 2) ) i E 9 2m(S_2) sin2 V mx + 

{m=0 



+^EE« 



(2c-d)(s-2) 



sing c xsinq c d xcos(q 2c - q 2( - c d) )t 



c=l d=l 



l-g^ s - z J ) ■( E T -''sin" </'"./■ + 






^E9 2c(s_2) sin^' 



y^g d(s 2 ^sinq c d a; cos 
d=l 



(M)' 9 2M E« 2 " 



Q = 






2m ( a " 2 > sin 2 j m i - 



rn— 



+- U-« 2 



2(s-2)\ V^ 2c(s-2) 



E* 2 



sin g a; • 



• E q- d(s ~ 2) smq c - d xcos{q 2 - l)^^ 1 ) + . . . + q 2 ^ d ^)t (56) 
d=l 

=:P x (x)+P xt (x,t). 
Note that the time-independent part 



P.M = i (l -««-*>) £ 



2m(s _ 2) l_cosg™2x 



m=0 



1 h _ «2(s-2)\ °° 

1 - iL_? 1 Y q^'~*) co S q m 2x, 

7T IT *• ' 



m=0 



is a Weierstrass-like function with the dimension s' = max{2s — 2,1} € [1,2) 
(i.e. for s 6 [1,3/2], s' = 1). From the equation (|5q ) one immediately gets the 
spectrum of P(x,t): all the frequencies governing the time evolution are 

u; c , d =(q 2 -l)(q 2 ^ + ...+ q ^- d y) 7 

where c=l,2,..., rf=l,2,...,c. Thus all the frequencies divide by q 2 — 1 which 
is also the smallest frequency, so the fundamental period of P(x, t) is 2ir/(q 2 — 1). 



A. 2. Average velocity. Let us study the behavior of (x) 



r = / dxx\\P[ 
>o 



2 7T l ' 



" ^Et^-^cos^-I)*. 



fe=i 



(q 2k - If 
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The above expression is valid only for even q. For odd q we have just the first 
term, which is tt/2. 

The average x(t) is of class C 1 , because its derivative is given by an absolutely 
convergent series: 



d(x) 



dt 



16 



fe(s-l) 



7( 1 -«*-' , )£!b=i-<«"-i> 



fc=l 



<cj> 



Jfe(s-l) 



-,2k 



fc=l 



,s-3 



2c- 



1 



(57) 



We showed in section (|3|) that (|57]) is fractal, while for odd q the average velocity 
\d(x)/dt\, of course, is not. This seemingly strange behavior is caused by the fact 
that 

dx sin nx sin mx 

is non-zero only for m, n of different parity. However, if one slightly disturbs 
our function, for instance by changing an arbitrary number of terms to the next 
higher or lower eigenstates, the dimensions D x and D t will not be altered, but 
the average velocity will become fractal. In other words, with probability one, 
independently of the parity of q, the average velocity of the wave function 



$ {x,t)=M 



i(s-i) 



n=\ 



sin((g n ± l)x)e 



-i(g"±l) 2 t 



(58) 



is fractal characterized by the same dimensions D x and D t as the function cur- 
rently studied. 

An explicit example of a similar function for odd q it is 



$ 1 (x,t)=M 1 



2 s " 2 sin(2a;)e" 



-i2*t 



Y" q n(s-2) 
ra=l 



sm(q n x)e 



-iq n t 



(59) 



One can see the only difference between this example and the original one (@) 
is in the first term. This difference accounts for the smoothness or roughness of 
the average velocity. It is very interesting because normally one expects that it 
is the asymptotic behavior that determines the fractal dimension. Here we have 
exactly opposite case: a change in the first term (varying most slowly) of a series 
changes the dimension of a complicated function (v) . 

Average velocity of the wave packet (59) is smooth for even q and fractal for 
odd q. A function which gives fractal average velocity for both even and odd q 
is 



<P 2 (x,t) = M 2 
= M 2 



! sin(2x)e 



oo 
n=0 



l(s -^ sin(q n x)e- 



■iq zn t 



2 S ~ 2 sin(2a;)e 



-%i z t 



N 



${x,t) 
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where M^ is the normalization constant. On the other hand, 



$ 3 {x,t) =M 3 



J2 1 n{S ~ 2) sm(q n x)e 



gives smooth average velocity for both even and odd q. 
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